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ABSTRACT 

Upcoming space-based photometric satellites offer the possibility of detecting continuum flux variabil- 
ity at the micro-magnitude (/imag) level. We show that the Doppler flux variability induced by the 
reflex motion of stars due to planetary companions has an amplitude of (3 — a)K/c, where K is the 
reflex radial velocity amplitude and a w d In F v /d In v is the logarithmic slope of the source spectral 
flux in the observed frequency band. For many of the known close-in planetary systems with periods 
P < 0.2 yr, the periodic Doppler variability, 0(/imag), is significant relative to the variability caused 
by reflected light from the planetary companion. For companions with P > 0.2 yr, the Doppler signal 
is larger than the reflected light signal. We show that the future photometric satellites should reach 
the sensitivity to detect this Doppler variability. In particular, the Kepler satellite should have the 
photon noise sensitivity to detect at a signal-to-noise ratio S/N > 5, all planets with minimum mass 
MpSini > 5Mj and P < 0.1 yr around the ~ 10 4 main-sequence stars with spectral types A through 
K and apparent magnitude V < 12 in its field-of-view. 
Subject headings: planetary systems, techniques: photometric 



1. introduction 

Planets outside the solar system orbiting main- 
sequence stars were first discovered in the mid-1990s 
(Mayor & Queloz 1995; Marcy & Butler 1996) through 
the small periodic reflex motion they imprint on the ra- 
dial velocity of their primary stars. The planetary sys- 
tems discovered by the radial velocity surveys to date 
are qualitatively different from the Solar System in im- 
portant ways that were unexpected by earlier theoretical 
models of planet formation. 

Although radial velocity studies have been highly effec- 
tive at identifying more than a hundred extrasolar plane- 
tary systems 3 , the method is inherently limited in several 
ways. First, it can generally only be used to find planets 
orbiting old, late-type stars. Second, it is not sensitive to 
planets with mass M p < M® . Finally, one can only infer 
the semi-major axis a, eccentricity, and a lower-limit to 
the planet mass. It is therefore important to seek new, 
complementary techniques that will allow better charac- 
terization of the physical properties of the known planets, 
and extend the current sample to different types of stars 
and lower-mass planets. 

The power of combining the radial velocity method 
with other techniques was best illustrated in the case of 
the planetary companion to HD209458. Originally dis- 
covered using radial velocity surveys, this planet was sub- 
sequently found to transit its parent star. Combination 
of the radial velocity data with the original ground-based 
transit data allowed the measurement of the mass (with- 
out an inclination angle ambiguity), radius and density 
of this planet (Charbonneau et al. 2000; Henry et al. 
2000). Further ultra-high precision HST/STIS observa- 
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tions of several transits allowed the exquisite measure- 
ment of these planetary properties, as well as the detec- 
tion of sodium lines in its atmosphere (Brown et al. 2001; 
Charbonneau et al. 2002). Perhaps the most exciting as- 
pect of the HST measurements was the demonstration of 
the ability to achieve photometric precision at the level 
of ~ 10~ 4 with space-based observations. 

In the future, several space-based photometric satellite 
missions are expected to reach a precision better than 
that demonstrated for HD209458, down to the micro- 
magnitude (/imag) level. These include Microvariability 
and Oscillations of STars (MOST) , COnvection, RO- 
tation and planetary Transits (COROT) 5 , Kepler 6 , and 
Eddington 7 . Scheduled for launch in June 2003, MOST 
is a targeted mission, and among its primary targets are 
five stars with known extrasolar planetary companions. 
CO ROT, Kepler, and Eddington, will monitor thousands 
of stars in one or several fields and search for very low- 
amplitude variability, as produced by stellar oscillations 
and planetary transits. 

At the /imag levels of variability probed by these 
satellites, many otherwise undetectable effects associated 
with planetary companions become significant. For tran- 
siting planets, the effects of the quadrupole moment of 
the parent star and more distant companions (Miralda- 
Escude 2002), the oblateness of the planet (Hui & Seager 
2000, Seager & Hui 2002), and the presence of moons and 
rings associated with the planet (Sartoretti & Schneider 
1999, Brown et al. 2001), are all potentially detectable. 
For both transiting and non-transiting companions, flux 
variations arising from tidal distortion of the parent star 
due to the planet (von Zeipel 1924, Drake 2003), back- 
heating of the star from the planet (Green et al. 2003), 
and reflected and scattered light from the planet's at- 

4 http://www.astro.ubc.ca/MOST/ 

5 http://smsc.cnes.fr/COROT/ 
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mosphere (Seager, Whitney, & Sasselov 2000; Sudarsky, 
Burrows, & Hubeny 2003; Green et al. 2003), may all be 
significant. 

Here we examine another potential source of periodic 
flux variability, namely the photometric flux variation of 
the primary star due to the Doppler effect. The reflex 
motion of the star induced by its planetary companion, 
causes slight Doppler variations in its photometric flux. 
In §2 we derive the amplitude of this effect. In §3 we 
compare the Doppler effect with other sources of stellar 
flux variability in the presence of planetary companions. 
In §4 we demonstrate that the Doppler variability should 
be detectable with the upcoming MOST and Kepler mis- 
sions. Finally, we discuss the systematic errors and in- 
trinsic variability in §5, and summarize our conclusions 
in §6. 

2. DOPPLER FLUX VARIATIONS 

A star moving at a radial non-relativistic velocity v r 
relative to the observer, obtains a Doppler shift to its 
bolomctric flux by an amount [see Eq. (4.97b) in Rybicki 
& Lightman 1979], 
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where c is the speed of light and the subscript denotes 
observed quantities in the absence of the source motion. 
If the emitted flux per unit frequency scales as a power- 
law in frequency, F„q oc v a , then the Doppler shift in 
the observed frequency v = vq( \ + v r /c) implies that the 
observer detects, 
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The Doppler effect vanishes for a = 3 because F„/v 3 is 
a relativistic invariant (corresponding to the phase space 
density of photons). Therefore, F v = F v q for a = 3 
irrespective of the source motion. Any other value of a 
leads to a net Doppler effect due to the change in the 
photon number flux. 

For a general source spectrum one may substitute a rs 
{dhxFvo/dlni'o), where the derivative is averaged over 
the observed band of frequencies. Assuming a blackbody 
source spectrum with a temperature T e g, 
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where x = hv/kT e g. In the Rayleigh-Jeans part of the 
spectrum {x -C 1), a = 2, whereas in the Wien tail 
(x >■ 1), a = 3 — x. For solar-type stars (T e g ~ 5700 K), 
and observations in the optical band (A w 600 nm or 
v w 5 x I0 14 Hz), we get a » —1.3. Given that the 
amplitude of the Doppler signal and the amount of flux 
received are both functions of frequency, the signal-to- 
noise ratio will also depend on the observed frequency. 
For Poisson errors, the noise variance will be proportional 
to the square-root of the flux. We find that the signal- 
to-noise ratio for T c g ~ 5700 K and fixed logarithmic 
frequency interval peaks in the V-band, at A ~ 500 nm 
or v ~ 6 x I0 14 Hz. 

3. APPLICATION TO PLANETARY COMPANIONS 

The reflex motion of the star due to a planetary com- 
panion should produce periodic variations in its observed 
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Fig. 1. — Expected flux variation amplitude as a function of 
the orbital period for the known extrasolar giant planets (filled 
dots). Also shown are the expected signals for fixed planet masses, 
assuming circular orbits, inclinations of 90°, and a solar-type pri- 
mary. The dashed line shows the level of expected reflected light 
variations for a Jupiter-radius planet, with Lambert sphere scat- 
tering (albedo = 2/3) and at full phase. The dashed-dot line shows 
the level of expected ellipsoidal light variations for a Jupiter-mass 
planet and solar-type primary. The solid vertical line is the Kepler 
mission lifetime (4 years). 



spectral flux. For circular orbits, the radial velocity of 
the star relative to the observer is simply v r = K sin 2tv4>, 
where is the phase of the orbit, and K is the velocity 
semi-amplitude, 



K = 28.4 m s" 
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1 yr 



-1/3 



M p sin i 

Mr 
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(4) 

Here M p is the mass of the planet, Mj is the mass of 
Jupiter, i is the inclination of the orbit, is the mass 
of the star, and we have assumed that M p <C M*. Thus, 
the flux variability signal has a period P and fractional 
amplitude AF/F = (3 — a)K/c. For orbits with an 
eccentricity e, the right-hand-side of equation Q should 
be multiplied by a factor (I — e 2 ) -1 / 2 , which is typically 
very close to unity. 

Since K is a direct observable for companions de- 
tected spectroscopically via radial velocity variations, we 
can predict the amplitude of the the flux variations for 
known planetary companions with essentially no ambi- 
guity. Figure I shows the expected flux variation ampli- 
tude as a function of the orbital period for all the known 
radial- velocity planets that have been detected so far, in- 
cluding the new transiting planet OGLE-TR-56 (Udal- 
ski ct al. 2002, Konacki et al. 2003). We have estimated 
T G fT values from the known spectral types of the host 
stars, and assumed observations in the V-band. Also 
shown are the expected signals as a function of period 
for fixed planet masses, assuming circular orbits, inclina- 
tion of 90°, solar-type primaries, and observations in the 
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F-band. 

The reflected light amplitude for a planet of radius R p 

is 



100 



AF = 



(5) 



Here p and R p are the geometric albedo and the radius 
of the planet, and a is the semi-major axis of its orbit. 
The dashed line shows AF/Fo for Lambert sphere scat- 
tering (p = 2/3), and planets at full phase. Averaged 
over all phase angles, Lambert sphere scattering repre- 
sents the upper limit to the amount of reflected light 
for a wide range of reasonable models, although larger 
scattered light signals are possible near zero phase angle 
(when the planet is at superior conjunction) due to, e.g. 
strong back-scattering (Seager et al. 2000). 

For short-period (P < 0.2 yr), massive (M p smi > 
Mj) planets, the Doppler variability is expected to be a 
significant contaminant to the reflected light signal. For 
known radial-velocity companions, the Doppler effect is 
completely predictable and its subtraction is straight- 
forward. When searching for planets via reflected light, 
the magnitude of the Doppler effect will not be known 
a priori. However, since the reflected light signal will 
generally be largest when the planet is at superior con- 
junction, it will be 90° out of phase with the Doppler 
signal, and it should be possible to cleanly disentangle 
the two signals. The Doppler signal is expected to dom- 
inate for long periods of P > 0.2 yr, since the reflected 
light signal falls off more rapidly with increasing orbital 
radius (oc a~ 2 ) than the Doppler signal (oc aT 1 / 2 ). 

There will also exist ellipsoidal light variations due to 
tidal effects on the star from the planet. We estimate the 
amplitude of these variations as, 



AF ^ Mp flU 
F ~ M* I a 
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where iZ* is the radius of the primary, and (3 ~ 0.45 
is the 'gravity-darkening' exponent for solar-type stars. 
The amplitude of the expected ellipsoidal variations as a 
function of period is shown in Figure f as the dot-dashed 
line for M p = Mj, M* = M , and R* = R Q . This effect 
should be subdominant for all but the shortest periods. 
Green et al. (2003) predict that back-heating of the star 
by the planet should be an even smaller effect, O(10" 7 ). 

4. DETECTABILITY 

For observations over a net duration T, the expected 
signal-to-noise ratio can be estimated as, 

S (3 — a) K I „ , . 

— = - — ^-vrx Tx io-°- 4V \ 7 

N V2 c 

where Y is the number of photons per unit time collected 
by the instrument for a star of magnitude V — 0, and 
we have assumed P <C T. Note that the \f2 factor is 
just the root-mean-square of the sinusoidal signal (for a 
planet in a circular orbit) in units of the amplitude. One 
can improve the signal-to-noise ratio by y/2 if the period 
and phase of the signal are known in advance by making 
observations when the signal is at maximum. 

The MOST satellite will collect Y ~ 1.6 x 10 8 photons 
per second for a V = star (Green et al. 2003). In Figure 
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Fig. 2. — Mass as a function of orbital period for the known 
cxtrasolar giant planets (filled dots). The large solid squares show 
those known planets for which the signal-to-noise ratio (S/N) of 
the flux variation is > 3 for 100 hours of integration time with the 
MOST satellite. The lines show the S/N > 5 detection limits for 
the Kepler mission for stars with V = 12, and spectral types B, A, 
F, G, and K (upper to lower diagonal lines). 



2, we show MpSini versus period of the known extraso- 
lar planets. Taking into account the ^-magnitudes of the 
parent stars, we have calculated the expected signal-to- 
noise ratio (S/N) for a total of 100 hours of integration 
with MOST. Since MOST is a targeted mission, we have 
assumed that the observations can be spaced over many 
periods and made at phases when the signal is largest. 
Thus we assume a signal-to-noise ratio that is a factor 
of -\/2 larger than given in equation 10. We have high- 
lighted those stars for which S/N > 3. Of these, one is 
currently included in the MOST target list: r Boo, for 
which S/N ~ 6 is expected. 

The Kepler satellite will look at a single field contin- 
uously for T — 4 years, and will collect Y ~ 1.4 x 10 10 
photons per second for a V = star. The detection limits 
for primaries of spectral type B-K are shown in Figure 2, 
assuming a total integration time of T = 4 years, V = 12, 
circular, edge-on orbits, and S/N > 5 for detection. 
There are ~ 10 4 such main-sequence stars in Kepler's 
field-of-view. If the frequency of massive (M > 5Mj), 
short-period (P < 0.1 yr) planets is 0.1%, K epler should 
detect ~ 10 such companions. 

Both MOST and Kepler will have broad wavelength 
filters in order to gather as many photons as possible. 
MOSTs bandpass will cover 350-700 nm, whereas Ke- 
pler's bandpass will cover 400-850 nm. The spectra of 
the sources will vary considerably over these wavelength 
ranges, and our assumption of a fixed spectral slope ap- 
propriate to the center of the 1^-band will break down. 
We have estimated the magnitude of this effect by inte- 
grating the blackbody spectrum over a uniform bandpass 
from 400-800 nm. We find that the expected signal-to- 



4 



PERIODIC DOPPLER FLUX VARIABILITY 



noise ratio is smaller by ~ 5% than we have estimated 
before. The source spectra will deviate from the black- 
body shape due to the presence of, for example, spec- 
tral lines. To determine the effects of such deviations, 
we have calculated the signal-to-noise ratio using a high- 
resolution Solar spectrum of Kurucz et al. (1984), kindly 
provided to us by E. Turner. Again, we find that the 
difference from the naive calculation is small, ~ 6%. 

5. SYSTEMATIC ERRORS AND INTRINSIC VARIABILITY 

The amplitudes of the Doppler effects we have consid- 
ered are extremely small, even for the extremely precise 
upcoming space-based photometric missions. We there- 
fore briefly consider both systematic errors and intrinsic 
variability. 

The photometric stability requirements for the detec- 
tion of the signals we consider are considerably more 
stringent than the requirements for the detection of the 
signals these satellites are designed to measure. To 
achieve the photon limit during the duration of the tran- 
sit of an Earth- like planet, the photometry must be stable 
to roughly 1 part in 10 5 over the ~ 10 hour duration of 
the transit. In contrast, we require that the photometry 
be stable to roughly 1 part in 10 6 at least over the ~ 3 day 
periods of the smallest planetary orbits. Reflected light 
signals also require stability over similar time scales, but 
in the cases where these signals are several times larger, 
the stability requirements are less severe. 

Not much is known about intrinsic variability at the 
levels we are considering for any stars except for the 
Sun. The intrinsic variability of the Sun is < 10 -5 on 
time scales of a few days or less (Batalha et al. 2001). 
This is of order the expected integrated Poisson error 
over these time scales, and thus should not significantly 
affect our signal-to-noise ratio estimates. However, the 
Sun does show significant, of order a few parts in 10 , 
variability on 30-day periods or longer (Batalha et al. 
2001). This variability is a consequence of magnetic ac- 
tivity and sunspots. Depending on the coherence of its 
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periodicity, this variability may complicate the detection 
of planets with longer periods via the Doppler variability. 
If measurements in two frequency bands are available, it 
may be possible to isolate the Doppler effect through its 
unique dependence on a (see Eq. |2J . 

6. CONCLUSIONS 

Equation J5J) provides the amplitude of the Doppler 
flux variability induced by the reflex motion of stars due 
to close-in planetary companions. For many of the known 
close-in planetary systems, the Doppler variability will be 
similar or larger than the variability due to reflected light 
from the planetary companion (see Fig. 1). 

We have found that future space-based photometry 
missions should reach the required photon-noise sensitiv- 
ity for the detection of this Doppler effect. In particular, 
MOST should detect this effect due to the known com- 
panion of t Boo at a signal-to-noise ratio of ~ 6 in 100 
hours of integration. Kepler should detect all planetary 
companions with masses > 5Mj and periods < 0.1 yr to 
the ~ 10 4 main-sequence stars with V < 12 in its field- 
of-view (see Fig. 2). This may uncover planets around 
massive primaries that are inaccessible to radial velocity 
surveys and too distant from their primary stars to be 
detected via reflected light. 
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